Long range superexchange - an exchange interaction through empty bands 
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We derive a generalization of the RKKY interaction to semiconductors using perturbation theory 
on a non-degenerated two-impurity Anderson model. In metals the interaction is mediated by 
excitations of free carriers over the Fermi-energy. In semiconductors, where no carriers are present, 
the only possible excitations are those of the localized impurity electrons (or holes) themselves. Thus 
a possible interaction is closely related to superexchange. We find an oscillating anti-ferromagnetic 
spin-spin coupling due to impurity electron (hole) excitations. By treating the coupling through 
empty bands (superexchange) along the same route as carrier mediated interactions (RKKY) it is 
easy to compare these two kinds of spin-spin coupling. The interaction derived here is of special 
interest for diluted magnetic semiconductors. 
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Usually superexchange is formulated within a cluster 
model consisting of three sites: two cation orbitals are 
partly filled thus forming an effective spin moment and 
one intermediate anion orbital is completely filled. In 
fourth order perturbation theory the resulting spin-spin 
interaction between the cation sites reads (180° Mn-O- 
Mn): 



J = 



2V 4 



AIcT 1 ), 



(1) 



where ATo is the difference between the ground state and 
a configuration where one electron is transferred from the 
anion to the cation. U is the on-site Coulomb interac- 
tion at the cation and V is the hybridization between 
both kinds of electrons. 

However it is out of question, that the superexchange has 
also a long range component. The latter is very impor- 
tant for diluted magnetic semiconductors, especially for 
doped (II, Mn)VI semiconductors. A competition be- 
tween RKKY iniej^aption and superexchange is typical for 
these materialsErau'la. Nearest neighbor superexchange 
leads to local spin singlets that reduce the effective con- 
centration of Mn spins. The superexchange between 
more distant pairs of Mn gives an anti- ferromagnetic cou- 
pling that competes with RKKY as soon as free carriers 
are present, which may be generated by doping with N 
in (ii, Mn)VI semiconductors or are present from the 
very beginning in Mn)V semiconductors. To get 

a qualitative picture of the interplay between the differ- 
ent exchange interactions it is convenient to have some 
simple limiting expressions of the mechanisms at hand. 
These expressions should give an idea of the dependence 
of the exchange mechanisms on certain model parame- 
ters like e.g. the intraatomic exchange coupling or the 
local Coulomb repulsion of Mn d-electrons and the (sp)- 
d hybridization. For superexchange Eq. ([!]) shows these 
dependencies. 

Concerning the long range part of the interaction ad- 
ditional informations are crucial, i.e. the dependence on 
the inter spin distance A and the dependence on the elec- 
tronic band structure of the host material. Both proper- 
ties can not be deduced from Eq. (|l|) or from expressions 
derived from any other cluster model. It is the intention 
of this paper to derive an expression for the long range 



component of superexchange, which is well comparable 
to the RKKY interaction. To this aim we will study a 
toy model and adjust the parameters in such a way, that 
the so to say "standard RKKY situation" is recovered. 
That means two isolated spins should be located at a 
certain distance in a host material that is described by a 
nondegenerate uncorrelated band. For this limiting case 
we will then apply fourth order degenerated perturbation 
theory. 

The paper proceeds in the following way: First the most 
important indirect exchange mechanisms are briefly re- 
viewed and discussed. We will concentrate on superex- 
change, RKKY, Bloembergen-Rowland interaction and 
a mechanism similar to the latter as well as to superex- 
change. Compact expressions for the last three coupling 
mechanisms are discussed. Then we will introduce the 
toy model and adjust the parameters in such a way, that 
we reach the best comparability to the RKKY expression. 
In the next step we will derive an expression for superex- 
change that will be exact in fourth order perturbation 
theory for the prepared model situation. Since the toy 
model establishes a well defined limit for more complex 
calculations, this expression can be used as a check for 
certain approximations. For [demonstration we will com- 
pare a work of Larson et.alE with our result. Further- 
more the result should give a vivid idea of the distance 
and band structure dependence of superexchange. To 
this purpose we evaluate the superexchange expression 
numerically for some simple model lattices. 



I. INDIRECT EXCHANGE MECHANISM 

Indirect exchange mechanisms, i.e. effective spin-spin 
couplings usually between local spins at cation sites me- 
diated by diamagnetic anionS|-where intensively discussed 
in the 50th. e.g. by AndersonQ, Goodenough and othersQ. 
The main goal of these studies was to understand mag- 
netism in insulators such as MnO and to justify the use 
of the Heisenberg model for this class of materials. The 
works where primarily concerned with the leading inter- 
action of spins in adjacent lattice cells and consequently 
a lot of cluster models where adopted. 
A different topic is the effective spin-spin interaction in 
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metals. Here the interaction is mediated by free carriers. 
In the language of perturbation theory these carriers are 
virtually excited over the Fermi energy, which results in 
a spin-spin coupling that oscillates in sign and altitude 
in dependence of the inter spin distance. Such an in- 
teraction is usually called RKKY coupling. It was, first 
proposed by Ruderman and Kittel for nuclear spinstj and 
later generalized to electronic spins. It is often discussed 
e.g. in heavy fermion systems. In contrast to indirect 
spin exchange in insulators RKKY interaction is formu- 
lated within a band picture and it's dependence on the 
spin-spin distance is well known. 

In semiconductors both mechanisms, i.e. virtual excita- 
tions of carriers and non-carriers, may be important and 
even compete with each other. However for a lot of ma- 
terials the restriction to spins of neighboring lattice cells, 
typical for insulators, is not a good approximation any 
more. Let us mention Europium chalcogenides, where 
at least the next nearest neighbor cell is important or di- 
luted magnetic semiconductors (DMS). Therefore a band 
formulation that is analogous to RKKY is also desired 
for interactions caused by excitations of non-carriers, like 
e.g. superexchange. 

Such interaction types are widely discussed especially 
for semiconductors. Let us start with Bloembergen 
Rowlands interaction, that is the band analogue to the 
process described for clusters in Ref.B in Eq.. (29;30). 
A valence electron is virtually excited at site 1, both 
the electron and the hole are transferred and recom- 
bine at site 2. The spin of the electron and the hole 
are coupled to local spins at site 1 and 2 by an in- 
traatomic inter-orbital potential. This interaction is be- 
lieved to be responsible for the magnetic interactions in 
Eu-chalcogenidealj and also discussed for DMSLl A sim- 
ilar interaction (Eq. 23 of RefE is discussed by Litvi- 
nov and Dugaev for (III,Mn)V DMsN. In the follow- 
ing we will call the this interaction "impurity induced 
Bloembergen-Rowland interaction" since the impurity 
electrons (Mn-d electrons in (III,Mn)V DMS) are vir- 
tually excited instead of valence electrons. 
For all couplings discussed so far there is a " standard ex- 
pression" usually derived in perturbation theory for some 
toy model that describes pure basic conditions for the re- 
spective interactions. 

For RKKY this set up consists of two spins which are 
locally coupled to an uncorrelated partly filled electron 
band by an intraatomic spin-spin interaction J p d- The 
same holds for impurity induced Bloembergen-Rowland 
interaction, just that the "spin" is now described by 
partly filled localized electron orbitals and the electron 
band is empty instead of partly filled. The incomplete 
filling may be due to a strong on-site Coulomb repul- 
sion. For the classical Bloembergen-Rowland interaction 
we need again two spins and a completely filled valence- 
as well as an empty conduction band. Again, the spin 
of the electrons are coupled to the local spins by an in- 
traatomic interaction. 

The resulting expressions for RKKY and Bloembergen- 



Rowland interaction in perturbation theory read (natural 

units): 

RKKY: 



J(A) = 



pd sr^ cos((fc' - fe)A) 



E 



(2) 



The sum runs over all k within the Fermi sphere and all 
k! that are located outside. For parabolic bands one finds 
for the Bloembergen Rowland interaction: 



J(A) 



^> 2 ATo 

7T 3 A 2 



K 2 (2r/r ) r = (2mAT )' 1 / 2 



(3) 

with the band gap ATo, the inter spin distance A 
and the reduced effective electron mass m. K 2 is 
the Mac-Donald function. (K 2 ~ 1/A 2 for r <C 
r , K 2 ~ A" 3 / 2 e 2A / r ° forr > r ). 
The same holds for the impurity induced Bloembergen- 
Rowland interaction, where ATo is now the energy 
difference between the impurity and the conduction 
band and m is now the effective electron mass in the 
conduction bandE3. In the next section we want to treat 
superexchange and derive a similar expression for this 
interaction. 



II. MODEL DESCRIPTION FOR 
SUPEREXCHANGE 



The model that can describe a "pure" version of the 
long range superexchange should be similar to the above 
mentioned models, especially to the model for RKKY 
interaction. Thus the competition between the latter 
and superexchange can be studied. Our model consists 
of two impurity sites with an effective spin moment. 
This moment is due to partly filled localized orbitals, 
which is realized by a strong on-site Coulomb repulsion 
U at the impurity orbitals. Furthermore there is an 
free electron band , described by the dispersion efe, that 
is energetically separated from the impurities by an 
energy ATo. The chemical potential is located between 
the impurity orbital and the band. The latter is thus 
completely empty in the unperturbed ground state. The 
impurity orbitals and the band are kinetically coupled 
by a local hybridization V. The latter will constitute 
the perturbation in the following calculation. This is a 
minimal set to study superexchange. Therefore other 
features, like e.g. intra-orbital Coulomb exchange, are 
not taken into account. 

The Hamiltonian for the described model reads: 



H — H o + H y 

i=l,2 



U 



i=l,2 



E e fe n L 
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FIG. 1: Schematic picture of conventional RKKY and virtual 
RKKY (superexchange) interaction for large on-site Coulomb 
repulsion U — > oo. 



H, 



V 



1=1,2 

E 



/l.C. 



(4) 



Let us note, that the construction operators can stand 
for holes or for electrons. If and create electrons the 
situation is closest to the usual interpretation of RKKY- 
interaction, i.e. that electrons are polarized by a local 
spin. For 5 = 1/2 the model (|) may describe RKKY as 
well as superexchange, if the band is partly filled instead 
of empty. In this case virtual excitations of band elec- 
trons over the Fermi energy contribute to RKKY, while 
virtual excitations of the impurity electrons lead to su- 
perexchange (see Fig. |l|). 

However in most cases superexchange is constituted by 
virtual excitations of holes instead of electrons (e.g. in 
the "classical" case of MnO). Thus in most cases the con- 
struction operators have to be interpreted as hole cre- 
ators and annihilators. To and efc are now energies for 
holes and U is the Coulomb repulsion between holes. The 
two interpretations of (|]) are connected via particle-hole 
transformation 



with the well-known results: 



Ml 



(5) 



Hence in the unperturbed ground state the impurity or- 
bitals arc filled with one electron (thus creating local 
spins). Furthermore the band is completely filled. Pos- 
sible excitations are from the band into the impurity or- 
bital with an excitation energy 



-T, 



Oh ' 



(8) 



This describes a situation where superexchange is exclu- 
sively mediated by filled (valence) bands. Such a situa- 
tion is not only found for the classical magnetic insulators 
(MnO), but also in (II,Mn)VI semiconductorsBB 
For (III, Mn)V semiconductors, too, the models of 
RKKY and superexchange are well comparable if the con- 
struction operators apply to holes, since the important 
carriers are holes in these systems. 

However the simple model (^) is quite general and does 
not only apply to DMS but to every situation, where vir- 
tual excitations of localized electrons are important. 
Let us now derive an expression for long range superex- 
change within this model. The most instructive way 
of considering the virtual processes leading to the spin- 
spin interaction is perturbation theory since one sums 
explicitely over all excited states. Treating the hybridiza- 
tion term as the perturbation we find that the free ground 
state is four-fold degenerated (with respect to the spin 
configuration). While calculating the energy corrections 
it is convenient to characterize the eigen-states of the 
free Hamiltonian Hq by their number of impurity elec- 
trons, which is a good quantum number of the free sys- 
tem. Further the following property of the perturbation 
Hy should be considered: 

If Hy works on a free ground state it changes the 
number of conduction (p— ) and impurity (d— ) elec- 
trons (holes) by one (while the total number of elec- 
trons (holes) is conserved). 

Due te-, this all odd energy corrections E a E a 3 ^ . . . 



fishfi. 



vamsntJ. In second order we find an energy contribu- 
tion, which does not affect the degeneracy of the ground 
state. 



TL = -(Ti e + U e 



l 0h 
f P 

u h 



= -e p 



ke 

U e = U 



(6) 



Let us discuss the following situation: The construction 
operators apply to holes and thus the hole energies T* h , 
e p h and \ih are fixed. Now let us consider the limit U — » 
oo. In this case the one-electron energy of the impurities 
Tq c goes to — oo while the energy of an doubly occupied 
impurity orbital stays finite (T$ e + U = -T$ h ). The same 
holds for the Bloch energies e p 



z ke . For superexchange we 
want to discuss the parameter constellation T$ h < /i/,. < 



e p kh and U 



oo(see Fig. hi). For electrons this means: 



TL « 4e < Me < (Tol + U) 



E^ 



V 2 



(9) 



The degeneracy is broken not until fourth order pertur- 
bation theory which gives an energy contribution E„\ 
Besides a constant term, which does not affect the spin 
orientation, this is given by: 



4 4) = E 



Ijb Tjc rjd zja 
tl r , ■ tli ■ n„ ■ a. a. j 



bed {^n ~ E m ){tj n 



(10) 



where if jf 



(E^\H V \E^). 



E^ ) is one ground state of 



the free system with the ground state energy En 



(0) 



2T$. 



(7) \E{b,c,d}) are excited eigenstates of the free system with 
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FIG. 2: Effective spin-spin coupling J as a function of the 
impurity distance A at different energy gaps. The distance 
is measured in units of the lattice constant. The band gaps 
(ATb — ¥f) range from 0.1 — 2.6 eV (from top to bottom). 
The impurities are located along [001] of a simple cubic tight 
binding lattice. Other parameters: band-with W = 6 eV, 
V=0.16W. Inset: J as a function of the ATb at (001), (002), 
(003). 



the energies E^ and Eo ■ The sum goes over all 
excited states. Due to the special shape of the perturba- 
tion potential Hy the latter states must have a certain 
number of excited electrons (holes) to get a nonzero en- 
ergy correction. Since we consider the limit U — > oo the 
excited electrons (holes) are in the conduction (valence) 
band. There is exactly one electron (hole) in the band 
in \E^) and \E^) and exactly two electrons (holes) are 

located within the band in \Ec°^). Thus |£Vj,,c,d}) can be 
expressed in terms of construction operators working on 
states with two impurity electrons \aa x /3a y ) (a, ft = 1 or 
2;a=Tor|). 



I^ 0) ) 



l^ 0) > 



Pl ll j 1 d i<?l\ i < j l aa 'x) 

(k20- 3 k 3 ai) 
pl 4<T7 d m a 8 \ma 8 /3a y ) 



d = (ki (77/3(78) 



(11) 



After tedious but straightforward calculations one ar- 
rives at: 



(Ei 0) \XYZ\Ei 0) )e^ 




FIG. 3: As FigJ^, but for a two-dimensional quadratic lattice 
with the band- width W = 4eV, V = 0.25 W. 



where 7 = w ^ Ni and 

X = Hvvt ai i 1 - n t 2 )Pp<yi H VPq rJ . i pta i 

Y = d ja5 (l - n^ 6 )eJt 5 

Z = Pta i Pra 3 H V p+ {l~n d ) Pva7 H V (12) 



In the sum the subscripts i ... m denote impurity sites (1 
or 2), while the subscripts o . . . v go over all lattice sites. 
The k— summations are over the first Brillioune-zone and 
(7i ... (78 are spin-subscripts. The phase factor (f> reads: 

$ = h(R - R P ) + h(Rq - Rr) + h{R s - Rt) + h(Ru - 

R v ). Performing the sum and introducing impurity-spin 
operators as usual 

st = 

s t +/ } = ^T/D^a/T) 

we finally can write the energy contribution in terms of 
an effective Hamiltonian of Heisenberg-form that works 



on the free ground state \GSa 



(0)\ 



= (GS^\H cS \GS^) with 
H eS = -J(A)Sx-S 2 . 



(13) 



We find for the exchange integrals J (A) (A is given in 
terms of the lattice constant): 



F(A) 



F(A) = 2cos((fc 2 - k 3 )A) +4cos((/ci - fc 2 )A) 

+ cos((fci - k 4 )A) + cos((fci + k 4 - k 2 - h)A) (14) 
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This effective spin coupling is of the anticipated order 
for small distances A. Due to the four-fold sum in 

we can not give an analytic expression for the asymp- 
totic behaviour of J(A). However, since excitations over 
the band gap are necessary we expect an exponential de- 
cay. In contrast to the "classical superexchange" where 
the particles fluctuate between the impurities and a single 
degenerated intermediate state, now the electrons may 
hop into different Bloch-states and still cause an effec- 
tive interaction. 

For the zero-bandwidth limit, i.e. e p k — Tq for all k, the 
k— sum in ( p~4| ) goes only over F(A). Since each cosine- 
function adds now to zero, the interaction vanishes in 
this limit. This is the correct result, because the sites 
are completely decoupled in the zero-bandwidth limit. 
The numerical evaluation of Eq. (|i~4|) gives always an anti- 
ferromagnetic interaction that declines with the distance 
and shows certain oscillations (see Fig. ^). The interac- 
tion gets even more important for systems with reduced 
dimensionality. For a two dimensional lattice the mag- 
nitude of the interaction increases approximately by a 
factor of five. This is seen in Fig. ||, where we used the 
same parameters for the nearest neighbor hoppings and 
the gaps between the band and the impurity-level as in 

Fig ¥n 

Eq. ( |14| ) gives an exact result in perturbation theory for 
a well-defined limit. Other treatments of long range su- 
perexchange that may involve more complicated models 
but also some additional approximations can be com- 
pared in the limit Uh — > oo with Eq. (|l4|). 
Let us demonstrate this with an example in literature 
that treats (77, Mn)VI semiconductors. In Ref.B Larson 
et.al. investigate electron and hole mediated superex- 
change and a special kind of Bloembergen- Rowland inter- 
action (negative local J). They apply a multi-band model 
with a realistic electronic structure, a local Coulomb re- 
pulsion U between Mn — 3d electrons and a hybridization 
between Mn ions and the host material. As explained in 
their paper, the five degenerated Mn — d orbitals can be 
modelled by a single orbital plus a factor that depends 
on the Mn ground state only. Thus only a single orbital 
is considered at each Mn site. The authors found that 
the superexchange caused by virtual excitations of two 



holes dominate. After applying the limits U — » oo, sin- 
gle nondegenerate valence band and local hybridization 
(V(k) = V) to their result, we want to compare them 
witruEq. The result of Larson et.al. (Eq.(4.4) of 

Ref.El) is written with electronic parameters. To com- 
pare it with our result we have to perform a particle hole 
transformation (|^) and apply the just mentioned limite 
and simplifications. Then the Mn-Mn exchange of Ref.B 
reduces to: 



= 2 



kk' 



V 4 cos(fc - k')A 
(T$-e k y{T$-e k ,) 



(15) 



This is quite close to the exact result in fourth order per- 
turbation theory (Eq.(p^). The remaining discrepancies 
seem to be a fair raice for the complexity of the model 
investigated in RefJj. 

Finally let us discuss qualitatively the influence of free 
carriers on superexchange and RKKY. If free carriers are 
doped into the band the virtual excitations of this car- 
riers over the Fermi energy lead to RKKY interaction. 
Since the energy gap is much smaller for these carriers 
the RKKY contribution should dominate in the sum (1C|). 
Furthermore, since the band is now partly occupied there 
are less virtual intermediate states for superexchange. 
This gives a vivid explanation for the fact that superex- 
change is suppressed by free carriers as e.fen worked out 
by Qimiao Si et.al. for CuO (Fig.l of RefB) 
In conclusion we have derived a simple expression for 
long range superexchange in a well defined toy model. 
This expression is useful for qualitative discussions and 
constitutes a limit, which can be used to evaluate approx- 
imations in more complex models. We have given an ex- 
ample of one such comparison for the case of (II, Mn)VI 
semiconductors, where the long range component of su- 
perexchange is very interesting. However, as in the case 
of RKKY, the physical picture developed here is quite 
general and is applicable to all problems where virtual 
excitations of two electrons or of two holes lead to an 
effective spin-spin coupling between these electrons or 
holes. 

This work was supported by the Deutsche Forschungsge- 
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